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We study the concurrence of entanglement between two quantum dots in contact to Majorana 
bound states on a floating superconducting island. The distance between the Majorana states, the 
charging energy of the island, and the average island charge are shown to be decisive parameters 
for the efficiency of entanglement generation. We find that long-range entanglement with basically 
distance-independent concurrence is possible over wide parameter regions, where the proposed setup 
realizes a “Majorana entanglement bridge”. We also study the time-dependent concurrence obtained 
after one of the tunnel couplings is suddenly switched on, which reveals the timescales for generating 
entanglement. Accurate analytical expressions for the concurrence are derived both for the static and 
the time-dependent case. Our results indicate that entanglement formation in interacting Majorana 
devices can be fully understood in terms of an interplay of elastic cotunneling (also referred to as 
“teleportation”) and crossed Andreev reflection processes. 
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I. INTRODUCTION 

Over the past few years, the prospect of realizing Majo¬ 
rana bound states (MBSs) in superconducting hybrid de¬ 
vices has attracted a lot of attention, see Refs. m for re¬ 
views. MBSs are predicted to emerge, for instance, as end 
states in proximitized one-dimensional Rashba nanowires 
in a magnetic Zeeman field, where signatures for their 
existence have already been reported in transport exper¬ 
iments [iHin]- For these topologically superconducting 
(TS) nanowires, MBSs are located at the respective ends 
of the wire, and the overlap between their wavefunctions 
becomes exponentially small for large nanowire length L. 
MBSs can then intuitively be thought of as spatially sepa¬ 
rated “half fermions”. They are under discussion as basic 
ingredients for “topological” qubits in quantum informa¬ 
tion applications, where information is nonlocally stored 
and, thus, should be quite robust against local decoher¬ 
ence processes. Systems with MBSs are also predicted 
to allow for non-Abelian Ising anyon braiding statistics. 
Since Ising anyons do not allow for universal quantum 
computation by themselves, it is essential to thoroughly 
understand the physics of Majorana states coupled to 
conventional qubits, where the latter can be realized in 
terms of nanoscale quantum dots [IIHIH]. 

Majorana bound states are described by self-adjoint 
operators, 7^ = yj, which anticommute with each other, 
{ijilk} = and with all other fermion operators. A 
pair of MBSs effectively forms a single nonlocal fermion 
with annihilation operator 

/ = (71+*72)/ v ^, (1) 

where the eigenvalue nj = 0,1 of the number operator 

fif = /V = *7172 + 1/2 ( 2 ) 


describes the state of the MBS pair [T]. For large L, the 
hybridization energy between the MBSs forming this pair 
becomes exponentially small, 

e/~exp(-L/0. (3) 

The length ^ refers to the spatial size of a MBS and is 
related to the superconducting coherence length in the 
proximitized nanowire [T]. For A —>■ 00, the MBS pair is 
then equivalent to a single zero-energy fermion. 

In this work, we consider the simplest case of a single 
MBS pair and study the concurrence of entanglement, 
C, between two single-level quantum dots m tunnel- 
coupled to the MBSs, see Fig. for a schematic illus¬ 
tration of the setup. Entanglement is a key concept in 
quantum mechanics and represents an essential resource 
for many quantum computation schemes m- Two spa¬ 
tially separated systems are said to be entangled when 
their quantum states cannot be described independently. 
One may naively expect that because each MBS amounts 
to half of the same fermion state, the MBSs themselves 
are strongly entangled. However, this is a meaningless 
statement since only the MBS pair has a well-defined 
state representation, see Eq. §• On the other hand, the 
nonlocality of the / fermion level suggests that the two 
dots coupled to the MBSs could be fully entangled. For 
a grounded superconducting island, where the charging 
energy Ec is negligible, this question has been studied in 
Refs. PTI|25) . and one finds that precisely the opposite 
statement holds true in the long-distance limit e/ —>■ 0, 
i.e., both dots are perfectly disentangled. This result 
is consistent with the absence of correlations among cur¬ 
rents flowing through normal-conducting leads in contact 
to different MBSs [TJ [26] . 

In what follows, we are especially interested in the pos¬ 
sibility of long-range entanglement generation, such that 
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Figure 1 . Schematic illustration of the setup. Large cir¬ 
cles correspond to two quantum dots (QDi and QD2), where 
gate voltages Vi,2 allow one to vary the energy levels £1,2. 
A Rashba nanowire is deposited on a superconducting island 
with charging energy Ec, where the proximity effect in com¬ 
bination with a Zeeman field induces a TS phase. The dimen¬ 
sionless gate parameter Ug (see main text) is proportional to 
a gate voltage Vg regulating the average charge on the TS is¬ 
land. Small (red) circles represent MBS operators 71 and 72, 
where the corresponding wavefunctions are located near the 
ends of the nanowire and the hybridization e/ corresponds to 
the distance L between both MBSs, see Eq. Tunnel cou¬ 
plings Ai and A2 connect the dots to the respective MBSs. In 
addition, the Josephson energy Ej describes the coupling to 
another grounded bulk superconductor (SC). 


the concurrence C{L) does not decay up to very long dot 
distance L |2T]. We shall determine the detailed condi¬ 
tions for the realization of such a “Majorana entangle¬ 
ment bridge” in the setup shown in Fig. In particular, 
we find that the charging energy Ec plays an important 
role in allowing for long-range entanglement generation, 
i.e., one should work with a mesoscopic floating (not 
grounded) superconducting island. We note in passing 
that standard notions of long-range entanglement refer to 
an exponential decay, C{L) ~ exp(—L/£), and merely re¬ 
quire that the lengthscale £ is long [IZIIIH]- In particular, 
for the topologically trivial limit of our setup, i.e., when 
the superconducting island does not host MBSs, such an 
exponential decay of C(L) is expected, see also Ref. [53] , 
where i is equivalent to the superconducting coherence 
length. The length-independence of C{L) predicted for 
the device in Fig. holds as long as the charging en¬ 
ergy, which scales approximately as Ec ~ 1/T in the TS 
nanowire geometry, remains the dominant energy scale. 
In particular, Ec should exceed both the thermal energy 
fcsT and the MBS hybridization scale £/ in Eq. 

Indeed, charging effects can change the physics in a 
profound way by coupling MBSs to the dynamics of the 
Cooper pair condensate on the island. In Ref. |3D], Fu 
studied the transmission through a setup as in Fig. 
but with normal-conducting leads replacing the dots. 
For large Ec and when the gate voltage is tuned to a 
Coulomb blockade conductance peak, i.e., the dimension¬ 
less parameter Ug regulating the average island charge is 
half-integer, only two degenerate island charge states are 
accessible and the model effectively maps to a spinless 
resonant tunneling problem. In a slight abuse of termi¬ 
nology, the resulting nonlocal phase-coherent transfer of 


electrons from one lead to the other has been dubbed 
“teleportation” (TP) by Fu [30], see also Refs. 132] . 
However, it has been argued that entanglement may be 
achieved more efficiently under Coulomb blockade val¬ 
ley conditions, i.e., when stays close to an integer 
[33li35] . As we discuss in detail below, entanglement is 
then a consequence of elastic cotunneling m processes. 
In our setup, elastic cotunneling has a strongly nonlocal 
character due to the underlying MBS realization of the / 
fermion, and therefore this process has also been referred 
to as “teleportation” [331131], see Sec. Ill A below. When 
more than two MBSs are present on the superconducting 
island, a nonperturbative version of the TP mechanism 
triggers the so-called topological Hondo effect [361143] . 
However, in the present work we shall focus on just one 
MBS pair, where topological Hondo physics is absent. 

Besides TP processes, however, it is well-known that 
entanglement can also be generated through crossed An¬ 
dreev reflection (CAR) processes. So far, CAR-induced 
entanglement in a system as shown in Fig. [^ has only 
been discussed for the grounded device, where Ec = 0. 
As mentioned above, when also e/ = 0, the dots will be 
perfectly disentangled, but for finite e/ (and hence not 
too large L), CAR processes generate entanglement both 
for the case of normal-conducting leads [2TH231 [26] and 
for the dot case at hand [Ml Hi] . This entanglement can 
be probed through the violation of Bell inequalities when 
several MBS pairs are present mi- In general terms, 
CAR refers to the splitting of a Cooper pair on the island 
which produces entangled electrons on different dots (or 
leads) [331133]. Similarly, in the reverse process, a Cooper 
pair is created by the combined in-tunneling of electrons 
from different sides. For normal-conducting leads, local 
Andreev reflection competes with CAR, since both elec¬ 
trons can be created (or annihilated) also in the same 
lead, and hence perfect entanglement is not possible [25] ■ 
For the single-level dots studied here, however, local An¬ 
dreev reflection is absent and perfect entanglement due 
to CAR processes may be realized. Below we address 
the fate of CAR-generated entanglement as well as the 
interplay of CAR and TP processes when interactions 
{Ec ^ 0) are important. 

In addition, we have also studied the time-dependent 
concurrence, C{t), found after one tunnel coupling is sud- 
dently switched on. Such a “quench” can be experimen¬ 
tally implemented by changing the voltage on a finger 
gate in close proximity to the respective tunnel contact, 
cf. Ref. [4]. By monitoring the time-dependent concur¬ 
rence, one can achieve a better understanding of the 
timescales on which entanglement is built up in such 
a setting [M]. For the setup in Fig. with Ec = 0, 
this problem has recently been studied in Ref. [25] • We 
here confirm and analytically explain the observations of 
Ref. [3S] for the noninteracting case. More importantly, 
we also address the entanglement dynamics for finite 
charging energy, and reveal the underlying timescales 
governing this interacting case. 

Before going into a detailed discussion, let us briefly 
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summarize the main results of our work: (i) We consider 
the full crossover behavior of the concurrence all the way 
from vanishing to large charging energy, (ii) allow for ar¬ 
bitrary rig, (iii) carefully address issues related to fermion 
parity, (iv) include the effects of a Josephson coupling 
Ej to an additional bulk s-wave superconductor, (v) ad¬ 
dress the L-dependence of the concurrence by including 
the direct overlap e/ between the MBSs, and (vi) study 
the timescales of entanglement generation after the sud¬ 
den change of a tunnel coupling. Throughout the paper, 
we provide closed analytical expressions for the concur¬ 
rence valid in different parameter regimes, both for the 
static and the time-dependent case, which are shown to 
give highly accurate approximations for our numerically 
exact results. Importantly, our results can be quanti¬ 
tatively interpreted in terms of an interplay of TP and 
CAR processes, thereby providing a comprehensive phys¬ 
ical framework to describe entanglement in MBS devices. 
In contrast to Ref. [33], we provide a rather complete 
analysis of this problem, including both TP and CAR 
processes, and providing detailed analytical results for 
the concurrence. 

The structure of this paper is as follows. In Sec. |ll 
we briefly summarize the model for the setup in Fig.^ 
For the detailed derivation of the model, we refer the 
reader to Refs. |3DH33] . We then turn to a discussion 
of the entanglement-generating tunneling mechanisms in 
Sec. |III[ The concurrence, C, provides a quantitative 
measure for entanglement of the two dots. We compute 
C analytically in several parameter regimes of interest, 
and validate the results by comparing to exact numerical 
diagonalization in Sec. |IV[ The time-dependent concur¬ 
rence obtained after a quench of tunnel couplings will 
then be studied in Sec. Vj Finally, we offer some con¬ 
cluding remarks in Sec. Wj Throughout the paper, we 
focus on the most interesting zero-temperature limit and 
use units with h = 1. 


II. MODEL AND DEFINITIONS 


A. Low-energy model 


We start by briefly summarizing the low-energy Hamil¬ 
tonian, H = Hyj + Hj + Hd + Ht, describing the setup 
shown in Fig.[^ see also Refs. [301133] . H contains (i) the 
part Hiij for the central TS island, where a mesoscopic su¬ 
perconductor with charging energy Eq causes proximity- 
induced pairing in the nanowire, (ii) the Josephson term 
Elj coupling the TS island to a second bulk superconduc¬ 
tor, (iii) the dot Hamiltonian El^, and (iv) a tunneling 
Hamiltonian Elt connecting the dots and the island. We 
consider energy scales well below the proximity-induced 
superconducting gap, such that it is justified to neglect 
quasiparticle excitations of the TS (see also the discus¬ 
sion in Sec. VI). The state of the island is then fully 


described by specifying the integer Cooper pair number 
Vc on the central island, where the number operator 


is canonically conjugate to the superconductor’s phase, 
[Nc,ip\- = —i, and the occupation number Uf for the / 
fermion built from the MBS pair, see Eq. 0. Including 
the MBS hybridization (§ and a Coulomb charging term, 
with the single-electron charging energy Ec and a dimen¬ 
sionless gate-voltage parameter Ug, the island Hamilto¬ 
nian is given by [301133] 

= Cf {hf - 1/2) + Ec (2Nc + nf - Ug^ , (4) 

where hf = f^f. The nonlocal / fermion representing 
the MBS pair thus couples to the condensate dynam¬ 
ics through Coulomb charging effects. Furthermore, the 
Josephson coupling Ej to another bulk superconductor 
(which is held at constant phase (po = 0 and different 
from the superconductor responsible for the proximity- 
induced pairing in the nanowire) is described by 

ffj = -Ejcos<p = -^(E‘^ + e-^‘^), ( 5 ) 

where the operator (e~*‘^) raises (lowers) the Cooper 
pair number by one unit, W —>■ IVc i 1. 

Next, each quantum dot in Fig. [^ is assumed to be in 
the Coulomb blockade regime such that the two dots can 
be modeled as single fermion levels at energies ei and 
£2, see Ref. [33]. These levels can be tuned by electro¬ 
static gates (Vj=i ,2 in Fig. and by the Zeeman field 
which induces the TS phase of the wire and also breaks 
spin degeneracy on the dots. As a consequence, we can 
use the effectively spinless fermion annihilation operator, 
dj, which is connected by a tunnel coupling Xj to 7j, see 
Eq.0 below. This tunnel coupling also captures a possi¬ 
ble spin dependence of microscopic transition amplitudes 
and can be taken as real-valued positive, see Ref. m- 
The Hamiltonian describing both dots then reads 

- 1 / 2 ), ( 6 ) 

i=i.2 

where the number operator fij = has eigenvalues 
rij = 0 , 1 . Equation ([^ represents a pair of charge 
qubits that can be entangled through the Majorana is¬ 
land. While charge qubits are more susceptible to detri¬ 
mental noise than spin qubits m. it is more transparent 
to analyze the concurrence within the present formula¬ 
tion. The extension of the above model to the case of 
spin qubits, where local qubit operations are also easier 
to implement, is rather straightforward by following the 
route sketched in Ref. [ 16 ]. However, this case is consid¬ 
erably more involved on a technical level, since one then 
effectively needs two TS nanowires. We here discuss the 
simpler spinless case, but due to the spatial separation 
of the dots, the results reported below represent “useful 
entanglement” [30] that could be exploited in quantum 
computation schemes. 

Finally, we come to the tunneling Hamiltonian iJj. 
Note that within our low-energy model, no TS quasi¬ 
particles are available for single-electron tunneling pro¬ 
cesses, and tunneling therefore has to involve MBSs. 
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We consider sufficiently large L such that di is tunnel- 
coupled only to the MBS described by 71 = {f + f^)/\/2, 
see Eq. 0” and Fig. Similarly, the MBS coupled to 
d 2 corresponds to 72 It is now crucial to 

take into account charge conservation, since a floating 
mesoscopic superconductor cannot simply absorb or emit 
charge 2 e without energy cost — this is possible only in 
the grounded case where Ec = 0. This consideration im¬ 
plies that tunneling terms ~ d'' must include a factor 
which annihilates a Cooper pair and thus restores 
charge balance. With the gauge choice in Ref. [3T], we 
obtain 

= [/ + + H.C., (7) 

V ^ i=1.2 

where “normal” tunneling terms like d^ f do not affect 
the condensate but “anomalous” tunneling terms such 
as d^e ~'^‘^change the Cooper pair number by one unit. 

In what follows, it will be convenient to choose occu¬ 
pation number basis states, 

|nin 2 n/,iV,) = (dj)"' (4)"' (4)”^ |000,iV,), ( 8 ) 

to represent the Hamiltonian. For this basis choice, the 
quantum numbers take the values rii^ 2 j = 0,1 and Nc = 
—-^max, ■ • ■, .^max, where fVmax restricts the number of 
Cooper pairs. Here A^max is taken relative to ng/2, and 
the A^max —> 00 limit of interest is rapidly approached 
when Ec is finite. The basis size is then 2M = 8(2A^inax+ 
1), and H becomes a 2M x 2M matrix. Diagonalization 
of this matrix yields the ground state I'k), which in turn 
determines the amount of entanglement between the two 
dots as explained below. 


B. Concurrence 

In order to quantify entanglement, one has to identify a 
suitable measure. For our bipartite case with two single- 
level quantum dots, a convenient and reliable measure is 
given by the concurrence, even though its construction 
is rather formal and abstract |20]. The concurrence is 
expressed in terms of the reduced density matrix for 
the dots, which follows from the full density matrix after 
tracing over the island degrees of freedom. In the zero 
temperature limit, this yields the Hermitian 4x4 matrix 

Pd = Tr„,,Arj4')(4'|, (9) 


where pj denotes the complex conjugate of pd and the Uy 
Pauli matrices act in the respective \nj) (with rij = 0,1) 
space corresponding to each of the two dots. Equation 
( 10 ) implies 0 < C < 1 , with C = 0 for separable (non- 
entangled) states and (7 = 1 for maximally entangled 
states. We note that the concurrence defined in Eq. (10) 


is a gauge-invariant quantity, i.e., we can use the above 
gauge choice with real-valued Aj in Eq. ([^. 

We have also computed the negativity m, which is 
an alternative measure for entanglement. While detailed 
results for the negativity differ from those for the con¬ 
currence, the physical conclusions are the same. From 
now on, we therefore discuss the concurrence only. Nu¬ 
merically, C follows from Eq. (10) once the full ground 


state of the total system, |4'), has been determined by 
diagonalization of the 2M x 2M Hamiltonian matrix. 

Analytical progress is possible when |4') has a simple 
form. An important example is given by 

1^) =Ci|000,A^l)+C2|011,Af2) + C3|101,Ar3)+C4|110,Af4), 

( 12 ) 

with arbitrary Cooper pair numbers Nj and complex co¬ 
efficients Cj subject to normalization, kiP = 1- 


Calculating C in Eq. (10) for this state, we find that 


at least one pair of number states in Eq. (12) must 


have identical (nj,A'c) in order to produce entangle¬ 
ment. There are two such possibilities, (i) Ni = N 4 with 
|ciC4| > IC 2 C 3 I, or (ii) N 2 = N 3 with IC 2 C 3 I > |ciC4|. In 
both cases, the concurrence is given by 


(7 = 2||ciC4|-IC 2 C 3 II, (13) 

see also Ref. [35]. Suppressing the {nf,Nc) indices, max¬ 
imally entangled states with (7=1 thus correspond to 
the Bell states 

|4>^) = ^(|01)+e“^-|10)), (14) 

|<1>b) = ^(|00)+e“^-|ll)), 

with arbitrary phases tpA/B- 


C. Parity conservation and symmetry relations 

Let us now address symmetries and conserved quanti¬ 
ties for our model. We first note that the total fermion 
parity defined as 


with trace equal to unity. The concurrence is then defined 
by HO] 

(7 = max 1^0, a/a)' - a/A^ - "vAs - vA) , (10) 

with the eigenvalues Ai > A 2 > A 3 > A 4 of the matrix 

G = Pd{(Ty®(Jy)pl{ay®CTy), (ll) 


is a conserved quantity since contains no terms mix¬ 
ing states with even and odd total electron number. The 
even-parity (odd-parity) sector has V = -1-1 {V = —1), 
where the parity-resolved ground states, I'k, P), follow by 
diagonalization of a Hamiltonian matrix of size M x M 
only. Conservation of V will generally be weakly broken 
in concrete physical realizations due to quasiparticles ne¬ 
glected in our model. In the presence of parity relaxation. 
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the true ground state then corresponds to the \ ^ ,V) state 
with lower energy. However, we here assume that V is 
conserved on all timescales of interest, such that both 
states ,V) are experimentally relevant and the con¬ 
currence depends on total fermion parity, C = C-p. Al¬ 
though this assumption represents an experimental chal¬ 
lenge due to the inevitable presence of residual quasi¬ 
particle poisoning, very recent experiments indicate that 
the corresponding timescales may reach minutes in simi¬ 
lar devices as considered here [45l [46] . 

We now turn to the n^-dependence of the parity- 
constrained concurrence, Cp(ng). Since a shift Ug —>■ 
Ug + 2 can be absorbed by shifting Nc —>■ Nc + 1 in the 
Hamiltonian, see Eq. Q, all observables have to be pe¬ 
riodic in Ug with period Arig = 2. (In the presence of 
parity relaxation, the period will in general be reduced 
to Aug = 1.) We thus obtain a first symmetry relation 
for the concurrence, 

Cv{ng)=Cv{ng + 2). (16) 

This relation allows us to restrict Ug to the window 0 < 
Ug <2 throughout. 

A second relation follows from an electron-hole-like 
symmetry property of H which relates the odd- and even- 
parity concurrences. 


Cp{ng]ei,e 2 ,ef) = (7-^(1 - ng;-ei,-£ 2 ,-e/). (17) 

In order to derive Eq. ( |l7| , we first note that H is 
invariant under the replacement £ 1 , 2 ./ with 

Ug —> 1 — rig, accompanied by a “particle-hole transfor¬ 
mation” U that exchanges creation and annihilation op¬ 
erators, / ^ etc., and —>■ yielding W — 

This symmetry implies that if we have the ground state 
|4^(ei,2./,’T'g))) we obtain another ground state from the 
relation 

|^'(-ei, 2 ./,l-ng)) = t7|4'(£i,2,/,ng)). (18) 

Now the point is that these two ground states have 
different total parities. Indeed, the original state has 

V = (—l)"i+" 2 -i-»i/^ while the transformed one has 

V = (_l)(i-’^i)+(i-"2)-£(i-»/) ^ _-p arrive 

at Eq. 0- In passing we note that for an island host¬ 
ing more than one pair of MBSs, if an odd number of 
dots is tunnel-coupled to MBSs, we instead would have 
V' = V, and Eq. ( [Tt] ) does not connect different parity 
sectors anymore. Moreover, the above arguments also 
show that the relative sign between and u/ in Eq. 0 
is irrelevant. 

For our island with a single MBS pair, Eq. 0 states 
that the concurrence in the odd-parity sector directly fol¬ 
lows from the even-parity result by simply inverting the 
energy scales £ 1 , 2 ,/ and letting rig —)■ 1 — Ug. From now 
on, unless noted otherwise, we therefore describe results 
for the even-parity sector V = +1 only, where it is suffi¬ 
cient to consider gate-voltage parameters in the window 

0 < rig < 2. 


III. ENTANGLEMENT-GENERATING 
PROCESSES 


Before turning to detailed results for the concur¬ 
rence of the two dots in Fig. let us first discuss the 
entanglement-generating processes at work in such a sys¬ 
tem, where the nonlocality of the / fermion built from 
the MBS pair turns out to be crucial. When searching 
for maximally entangled states of the two dots, Eq. (14) 


suggests two candidate sets for suitable superpositions of 
dot states \n 1 n 2 ), 

A = {|10),|01)}, H = {|00),|11)}. (19) 

In this section, we shall focus on the large-Ec case, as¬ 
suming that rig is not half-integer such that the equilib¬ 
rium island state is uniquely defined. Hence the Cooper 
pair number corresponds to A^ = 0, and for V = -1-1, we 
have n/ = 1 (n/ = 0) for set A (B). 


To realize the maximally entangled Bell states (141, 


(virtual) tunneling processes through the island now have 
to provide the necessary coupling between states within 
a given set. The resulting superpositions should ideally 
have equal weight, which is possible when the dot en¬ 
ergy levels are adjusted to fulfill the condition £1 = £2 
(£1 = —£ 2 ) for states in set A (B). Perturbation theory in 
the tunneling amplitudes then shows the existence of two 
different entanglement-generating mechanisms, namely 
elastic cotunneling (i.e., TP) and CAR. In fact, our quan¬ 
titative results will find a natural interpretation in terms 
of these two mechanisms. 


A. Teleportation 

The TP mechanism refers to the nonlocal transfer of 
an electron across the TS island. For example, let us 
consider the state |10), where an electron resides in the 
left dot and the right dot is empty. We now study how 
a coupling to the partner state |01) in set A is estab¬ 
lished by Ht- After the in-tunneling process from the 
left dot, the island charge will change by one unit. How¬ 
ever, this change typically comes with an energy cost of 
order Ec and, therefore, is possible only as a virtual pro¬ 
cess. The corresponding out-tunneling event restores the 
equilibrium charge state again and thereby may trans¬ 
fer an electron to the right dot. In effect, this process 
is similar to elastic cotunneling m but with two major 
differences. First, TP has a highly nonlocal character in¬ 
herited from the nonlocality of the / fermion. Second, 
because of the existence of “normal” and “anomalous” 
tunneling terms in Ht, there are two different contribu¬ 
tions. Using the notation in Eq. the first one (de¬ 
noted by TPji) proceeds solely by “normal” tunneling 
(^ d^f,pd), while the second contribution (TPq) only 
employs “anomalous” tunneling processes, 

TP„ : |101,0) ^ |110,0) ^ |0I1,0), (20) 

TPa ■■ |101,0) ^ |000,I) ^ |0I1,0). 
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By virtue of the processes in Eq. (20), the |10) and |01) 


states are coupled together. For ei = £ 2 , both states en¬ 
ter the superposition with equal weight. They may then 
form the maximally entangled Bell state Id)^) in Eq. (14), 
where TP is responsible for long-range entanglement. 

Finally, we note that the TP phenomenon discussed 
by Fu |30j is a resonant version of the above process, 
which takes place near half-integer values of Ug. We will 
address this case in detail in Sec. IV A[ where we show 
that perfect entanglement is not possible under Coulomb 
blockade peak conditions. 


B. Crossed Andreev reflection 

CAR provides a distinct second entanglement¬ 
generating mechanism, where for ei = —e 2 , the two dot 
states in set B enter the superposition with equal weight 
and may yield the Bell state |$b) in Eq. (14). Since en¬ 


tangled states have to involve a pair of number states 
with matching {nf,Nc) entries, see Eq. ( |I^ , the two 
states of interest now correspond to |000,0) and 1110,0). 
Their total particle numbers are different, and hence not 
only the tunnel couplings Ai_ 2 ) but also a finite Josephson 
coupling Ej to another bulk superconductor is needed in 
order to yield non-zero concurrence. We now show that 
this condition arises even though CAR processes can con¬ 
nect states with different island particle number 2 Nc + nf 
already for Ej = 0. Indeed, starting from |000, 0), there 
are six different CAR transitions corresponding to all 
permutations of the three elementary transfer processes 
^ Ai, A 2 , Ej. As an example, we show two of these CAR 
sequences, 

CARi : |000,0) ^ |011,-1) ^ |011,0) ^ |110,0), 

CAR 2 : |000,0) ^ |000,1) ^ |101,0) ^ |110,0), 

( 21 ) 

where the first (second) step for the CARi (CAR 2 ) pro¬ 
cess describes anomalous tunneling, the respective third 
step corresponds to normal tunneling, and the second 
(first) step involves the Josephson coupling in Eq. ([^. 
In the latter step, a Cooper pair is transferred to the TS 
island. It is only because of this step that the two states 
|000, 0) and |110,0) can be connected. Finally, we note 
that in our quantitative analysis below, not only the two 
processes in Eq. (21) but all six CAR sequences will be 


taken into account. 


IV. MBS-MEDIATED ENTANGLEMENT 

In this section, we quantitatively discuss the concur¬ 
rence C, describing the amount of entanglement between 
the two dots in Fig. by comparing analytical expres¬ 
sions in different parameter regimes to the correspond¬ 
ing numerically exact results obtained from the ground 
state of the full system, see Eq. (10). Since the entire 


behavior of the concurrence in the odd-parity sector fol¬ 
lows from the even-parity results, see Eq. ( [T7| , we shall 
only discuss the case V = -1-1. Furthermore, because 


of Eq. ((16|), the gate parameter Ug can be taken in the 


window 0 < Ug <2. 


For convenience, we will employ a compact notation 
for the number states by writing Iv) = \nin 2 'nf, Nc) in 
Eq. Without the Josephson coupling and the tunnel 
amplitudes, i.e., for Ej = Ai ^2 = 0, the eigenenergies of 
the full system are then given by 


El, — Ec{Q Tig) + ej(^i 1/2), 


( 22 ) 




with the integer island particle number Q = 2Nc + rif. 
The dependence of Ei, on Q includes the well-known 
Coulomb charging energy parabola m and is sketched 
in Fig.[^ We shall return to Fig.j^below when discussing 
our analytical results for the concurrence. Moreover, the 
symbol JV appearing in some equations below denotes a 
normalization constant for the respective state. 

Before describing our results for the concurrence in de¬ 
tail, we pause to offer some guidance for focused readers. 
In Sec. |IV A[ we analyze the limit of strong Coulomb 
blockade, where the charging energy is the dominant en¬ 
ergy scale and only a few island charge states are acces¬ 
sible. We derive analytical expressions for the concur¬ 
rence C for gate-voltage parameters near the Coulomb 
blockade valleys centered around Ug = 1, see Eq. (30), 


and ng = 0 mod 2, see Eq. (|36|) below. In addition, we 


also provide closed results valid near a Coulomb blockade 


peak, see Eq. (42), and compare our analytical results to 
full numerical diagonalization. We then briefly discuss 


the noninteracting {Eq = 0) limit in Sec. IV B where 
the concurrence follows as given in E q. ([48 ). Finally, we 
return to the strong-Ac limit in Sec. |IV C but now also 
allow for a strong MBS hybridization e/, i.e., arbitrary 
length L in Eq. However, for simplicity, we only 


study the case of integer Ug in Sec. IV C 


A. Strong Coulomb blockade regime 

We start our quantitative analysis of the concurrence 
with the limit of large Eq, 


Ec » max|ei,2,/, Ai,2,Ej|, 


(23) 


where the charging energy dominates over all other en¬ 
ergy scales except for the TS gap, which implicitly 
sets the ultraviolet cutoff for our model. The param¬ 
eter regime in Eq. (23) has also been addressed in 
Refs. It is convenient to distinguish four differ¬ 

ent, partially overlapping, gate-voltage parameter win¬ 
dows defined by 


1 3 

= ng.o + ??/2, e 0, -, 1, - 


( 24 ) 
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Figure 2. Schematic illustration of the energy E^, see 
Eq. (221, vs the island particle number Q = 2Nc + n/ for 
ei _2 = 0. Panel (a) shows the case Ug = 1 for Ec > £/• Panel 
(b) is also for Ug = 1 but with Ec < £/, while panel (c) is for 
Ug = 0. The respective even-parity states \v} = |nin 2 n/, A^c) 
are also indicated, where tunneling (A 7 ^ 0) and/or Josephson 
(Ej 7 ^ 0) processes cause the shown transitions. 


where [yyl < 1 parametrizes the deviation from the re¬ 
spective center value Ug = Ug^. For both half-integer 
values of Ug = Ug^, the dominant charging energy con¬ 
tribution, given by the first term in Eq. (22), is precisely 
degenerate and thus corresponds to a Coulomb blockade 
conductance peak |32]. We present our results for this 
case later on, but first turn to the Coulomb blockade 
valley at ng_o = 1- 


Coulomb blockade valley near Ug = I 


For rig 0 = 1 in Eq. (24), we observe from Eig.j^a) that 
the ground state has island occupation number Q = 1 in 
order to minimize the large Coulomb energy, such that 
Uf = 1 and Nc = 0. For V = -|-1, one additional electron 
then has to occupy the dots. As a consequence, the low- 
energy sector is spanned by 


Ho = {|101,0),|011,0)}, 


(25) 


with the two states denoted by Ivq). The energetically 
closest set of states, \vi) G 'Hi, has Q = 0 or Q = 2, 
see Fig. [^a), where the /-electron level representing the 
MBS pair is empty, nj = 0, with = 0 or = 1. The 
two dots are then either both occupied or both empty, 
resulting in the subspace 


Hi = {|000,0), |110,0), |000,1), |110,1)}. 


(26) 


These states are separated from the Ho sector by an en¬ 
ergy of order Eq- This separation of energy scales allows 
us to derive a low-energy Hamiltonian, H, describing the 
parameter regime in Eq. (231 for gate-voltage parameters 
close to rig^o = 1- This reduced Hamiltonian is obtained 
by projecting the full Hamiltonian to the subspace Hq 
only. Using the two basis states [t'o) in Eq- (25), H cor¬ 


responds to a 2 X 2 matrix, such that it is straightforward 
to determine the ground state and the concurrence ana¬ 
lytically. 

The above projection is implemented by disregarding 
all high-energy states beyond \vi) G Hi in Eq. (26) and 


by treating the coupling between the subspaces Ho and 
Hi through a Schrieffer-Wolff transformation [Ml 
which here is equivalent to second-order perturbation 
theory. For Eq > £/, this coupling is due to tunnel¬ 
ing (HU only, see Fig.j^a), and with the energies E^ in 
Eq. d^, we obtain 


g (^ c 

1- (HI), 

h = J 2 e. 


Wo){vo\ + 


(27) 


^0 


E 




1 


E,,„ — E,, 


E,a - E,, 


{vo\Ht\vi){vi\Ht\v'^), 


where the second term describes virtual excursions to the 
energetically higher |i^i) states. Using Eq. (23), £ 1 , 2 ./- 


terms in the respective denominators are small against 
the large charging energy contribution, and up to an over¬ 
all energy shift, we arrive at 


H = 

(a b 
b —a 

a = 

ei - £2 

2 

h — 

A 1 A 2 


Ecil- 


(28) 


2Hc(l-r7')’ 


where r] parametrizes Ug around Ug^ = 1, see Eq. (24). 
The ground state of Eq. (28) is 


|4i)=A7 (a-|101,0)-f 6|011,0) 


(29) 

































which matches the general form in Eq. (12). The concur¬ 
rence then directly follows from Eq. (13), 


C{X,) = 


1 


\/^ 


-xl 




(30) 


stating that C is a universal function of a single param¬ 
eter Xi which depends on basically all the microscopic 
parameters according to Eq. ( [^ . 

Eor Xi = 0, Eq. (30|) yields the ideal value, (7 = 1, 


characterizing perfect entanglement. Indeed, Eq. (29) 


then reduces to the Bell state |$^) in Eq. (14), and we 
conclude that entanglement is here established by the 
TP mechanism. This condition for perfect entanglement 
is fulfilled for all symmetric systems, i.e., for identical dot 
level energies and tunneling strengths. 


Cl — £2 — e, Ai — A 2 — A. 


(31) 


The approximations leading to the reduced Hamiltonian 
(28), and hence to the concurrence in Eq. (30), eventually 
break down near a Coulomb peak with \rj\ —>■ 1. Nonethe¬ 
less, Eq. (30) shows that (7=1 can persist over an ex¬ 


tended gate-voltage parameter window around Ug^ = 1. 


Coulomb blockade valley near rig = 0 

In a similar manner, we next study what happens 
around ng_o = 0. In Fig. [^c), we show a subset of 
the relevant states, implementing two of the possible six 
CAR processes. The energy is then minimized for 
Uf = Nc = 0, and the V = +1 low-energy sector is 
spanned by 


.^K,o-0) ={|000,0),|110,0)}, 


(32) 


where both dot levels are either occupied or empty. The 
excited states relevant for CAR processes, cf. Sec. |IIIB[ 
then correspond to 


0 ) 


= {| 110 ,- 1 ),| 101 ,- 1 ),| 011 ,- 1 ), 

|101,0), 1011,0), 1000,1)}, (33) 


where the first and the last state have island charge Q = 
=F2 and thus are highest in energy, ~ iEc, while the 
remaining states (with (J = ±1) are separated from the 
Ho sector by an energy scale ~ Eq- 

The two states in Eq. ( [3^ have different total particle 
number. As discussed in Sec. Ill B see also Fig. [^c), the 
CAR mechanism connecting those states has to involve 
a finite Josephson coupling Ej and therefore constitutes 
a third-order process. In particular, now a transition 


within the high-energy sector in Eq. (331 is necessary. 


and we thus have to include one normal tunneling, one 
anomalous tunneling, and one Josephson process, i.e., we 
have to apply iJt twice and Efj once, see Eqs. ([^ and ([^. 
The corresponding third-order perturbation theory can 
again be implemented by a Schrieffer-Wolff transforma¬ 
tion SZl, which effectively yields -f 77*'^^ for the 


reduced Hamiltonian in Eq. (27). With Elt = Ht + Hj, 
we find |48] 


{vo\Ht\vi){vi\Ht\v'i){v[\Ht\vo) 

(Ei/g — E^-^){E^g — E^i^) 


Wo){’' 


ol> 


77(3) ^ ^ 

. 

where \i'i) and |z/() are the virtually occupied higher- 
energy states in Eq. (33) and terms ~ £ 1 , 2 ,/ are small 


against the charging energy contributions in the denom¬ 
inator. Using the two states in Eq. ( [^ to represent 77, 
we obtain the same 2x2 matrix as in Eq. (28) but with 
the replacements 


a —>■ a = — 


b-^b' = - 


£1 + £2 7 ?(Af -I- A2) 

2 2 Ec{l-v^y 

3 X 1 X 2 EJ 

£;2(i_^2)(4_^2)’ 


(35) 


where ry now parametrizes rig around rig^o = 0- We ob¬ 
serve that second-order processes only give corrections to 
the state energies (o'), while third-order CAR processes 
produce the crucial coupling of both states (5') needed 
for entanglement. We note in passing that second-order 
corrections ~ Ej only give an irrelevant constant energy 
shift that has been dropped in Eq. (35). The ground 


state finally follows as in Eq. (29), but with the ba¬ 
sis states in Eq. (1^ and the replacements in Eq. (|^|. 


Hence the concurrence is again given by Eq. (30) but 
with Xi — Xq = a'/b'. 

Although these results look very similar, there are cru¬ 
cial differences in the entanglement properties when com¬ 
pared to the rigfi = 1 case. In particular, although we 
may have Aq = 0 for a specific gate-voltage parameter 
(ry), such that perfect entanglement is realized, it is not 
possible to have (7 = 1 over a wide parameter range any¬ 
more. In that sense, within the large-Tlc' regime defined 
in Eq. (23), entanglement produced by the CAR mecha¬ 


nism is less stable than the one caused by TP. To demon¬ 
strate this point, consider dot energies £1 = —£2 where 
CAR processes are most effective, cf. Sec. |HIB[ keeping 
symmetric tunnel couplings as before, Ai = A 2 = A. We 
then obtain 


CiXo) = 


1 




^ (4-ry2)ry7;c 

Ao - (36) 


describing a narrow concurrence peak centered around 
rig = 0, with height (7(0) = 1 and width brig ^ Ej/Ec- 
This concurrence peak is a clear signature of CAR pro¬ 
cesses and can easily be distinguished from the broad and 
robust (7 = 1 plateau generated by TP around rig = 1. 
In particular, the CAR entanglement peak quickly dis¬ 
appears with increasing \ei + £ 2 ], while it remains sta¬ 
ble against variations with £1 « —£ 2 - We thus obtain 
precisely the opposite behavior as for the TP-generated 
concurrence plateau. 
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Concurrence near Coulomb bloekade peaks 

Next we study the concurrence behavior near half¬ 
integer values of rig. For simplicity, we discuss the most 
interesting symmetric case defined in Eq. (31). With rig 


near rig^ = 1/2, there are four low-energy states corre¬ 
sponding to the subspace 


% 


("s,0 = 1/2) 


= {|000,0),|110,0),|101,0),|011,0)}. (37) 


Since these states are already directly coupled by tunnel¬ 
ing, the reduced Hamiltonian H can simply be obtained 
by projecting H to the subspace in Eq. (37), thereby 


neglecting all virtual excursions to higher energy states. 
When represented in the basis (37), up to an irrelevant 


overall constant, H then takes the form 

/E-2e 0 0 0 \ 

fring,0 = 1/2) ^ 0 E \ly/2 -Xly/2 

0 A/V2 -E 0 

0 -X/V2 0 -E J 


with the rig-dependent energy scale 

E = Ec{ng - 1/2)+ {e-ef)/2. 


The ground state of H is either given by 

^^2 + x^-E^ | 110 , 0 ) 


=A7 


(38) 


(39) 


(40) 


-A|ioi,o) +A|oil,o) 

or it corresponds to the separable (7 = 0 state 11112 ) = 
|000, 0). The respective energies are 


El = -Ve^ + X^, 


Eo — E — 26. 


(41) 


SinceKki) is of the form in Eq. (12), we can again use 


On the other hand, it can happen that at some gate 
parameter rig ' (above no for e > 0), the turnover 
value (7 = 1/2 separating the CAR-dominated regime 
(near the rig = 0 valley) from the TP-dominated regime 
(near rig = 1) is reached. Using Eq. (42), this point is at 
E = 0, corresponding to 


(C=i/2) _ 1 e/ - e 

® ~ 2^ 2Ec ' 


(44) 


Note that this expression holds in particular for e < 0. 
We therefore see that the gate parameter Ug^c determin¬ 
ing the transition through (7 = 1/2 is given by 


Un c = max 


(noc 


,( 0 = 1 / 2 ) 


_ 1 £/ + |e 

2 2Ec 


(45) 


where the last expression becomes exact for A ^ |e|. 
For no > ng^~^^‘^\ the concurrence transition is an 
abrupt jump. Otherwise, it is a smooth transition, and 
the concurrence jump only happens at a gate parameter 
Ug = no < rig^”^^^^ where (7 has already dropped to a 
small value. 

The actual value for Ug^c in Eq. ( |4^ is close to Ug = 1 /2 
for large Ec, but it exhibits a systematic shift towards 
Ug = 1 when |e| increases. We also note that as Ug 
moves from rig^c towards Ug = 1, the concurrence given 


by Eq. (42) approaches the perfect entanglement regime, 
(7 —>■ 1, as predicted by our expansion around Ug = 1. 
Consequently, our results for C{ng) obtained within dif¬ 
ferent Ug windows smoothly match onto each other. 

Finally, a similar analysis can be carried out for rig 
near ng_o = 3/2. The result is given by Eq. (42) again 


but with the replacements rig —>■ 2 — iig and e —>■ —e. 


For A < |e|, the transition corresponding to Eq. (45) now 
takes place at the gate parameter 


n„ r — - 

g,c 2 




2 Ec ’ 


(46) 


Eq. ^ to directly read off the concurrence, but E 2 > E^ towards rig = 1 with increasing |e|. 


is necessary to have a finite concurrence. With the Heav¬ 
iside step function 0, we thus obtain the concurrence in 
analytical form. 


Discussion 


C{E/X,e/X) 


1-f 


e(E2-Ei) 

(\/E2 + A2 - e) /a 


(42) 


which now depends on two parameters. In order to obtain 
a non-zero value in Eq. (42), the gate parameter must be 
above a critical value, rig > rig, with 


1 e/ + e A^ 

2 2Ec ~ 4:eEc' 


(43) 


This relation is only relevant when e > 0, since the con¬ 
dition E 2 > El is always met for e < 0. For e > 0, the 
concurrence thus exhibits an abrupt jump from (7 = 0 to 
a finite value as rig increases through ng. 


With analytical results for the concurrence at our dis¬ 
posal, we now compare them to numerically exact results. 
The rig-dependence of the concurrence for e/ = Ej = 0 
is shown in Fig. where we consider a symmetric sys¬ 
tem with e = 0. For Ec = 20A, the analytical curve 
almost perfectly matches the numerically exact result, 
where we have (7=1 for practically the full valley range 
1/2 < rig < 3/2, with a sharp crossover around rig = 1/2 
to (7 = 0 near even rig. For Ec = 5A, the qualitative 
behavior is still captured by our analytical expressions 
but there are quantitative differences, as expected when 
moving away from the strong-Ec limit. We note in pass¬ 
ing that for odd parity, the concurrence behavior in both 
valleys will effectively be interchanged. 
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Figure 3. Concurrence C vs Ug for the even-parity sector Figure 5. Same as Fig. but for finite e, using Ej = 1. 

with Cf = Ej — 0. We consider a symmetric system, see 
Eq. ( |31| l, with e = 0 and A = 1, for two different values of 
Ec- Solid curves are obtained by numerical diagonalization of 
the full model. Dashed curves represent the analytical result 
derived in the large-Sc limit, see Eqs. (1301) and (42l. With 


C(ng) = Cijig -b 2), see Eq. (161, we only show one period. 


c 

1 

0.8 

0.6 

0.4 

0.2 

0 



0 


0.5 


1.5 


Ur. 


Figure 4. Same as Fig. but for finite Ej and Ec = 20. 
The dashed curves again give analytical results, using Eq. (36) 
near = 0 and Ug = 2. 


Let us next address the effects of adding a small but fi¬ 
nite Josephson coupling Ej, where we again study a sym¬ 
metric system with e = £/ =0. The concurrence is shown 
in Fig. I^for Ec = 20A. We first notice that near Ug = 1, 
the effects of Ej are negligible and the TP-induced per¬ 
fect value (7 = 1 is robust against adding a small Joseph¬ 
son coupling. However, we now encounter the predicted 
narrow concurrence peak with width Sug ~ Ej/Eq cen¬ 
tered at Ug = 0 mod 2, 


see Eq. (36), which provides a 


very characteristic signature of CAR processes. 

In Figs. and 1^ we have considered a symmetric sys¬ 
tem with e = e/ = 0, where the concurrence exhibits 
the additional symmetry C{ng) = C{2 — Ug), and the 
turnover between both concurrence regimes is perfectly 
smooth and happens precisely at half-integer Ug. How¬ 
ever, both of these features will not be present in general, 
and a counter-example with finite e is shown in Fig. 

In accordance with Eqs. (45) and (46), we find that the 


regimes move closer to = 1 with increasing e. As ex¬ 
pected from our above discussion, we also confirm that 
for e > 0, the crossover in C{ng) (at ng slightly above 
Ug = 1/2) deforms into an abrupt jump, while the corre¬ 
sponding crossover near Ug = 3j2 stays smooth. Hence 
Fig. [^clearly shows that C(ng) ^ C{2 — Ug). We men¬ 
tion in passing that for finite but very small e, the narrow 
entanglement peak seen in Fig. may survive, but with 
the center position shifted to finite Ug. However, this 
peak is found to quickly disappear with increasing |e|, 
cf. Sec. |HI B[ and is altogether absent for the parameters 
in Fig. |5j 


B. Noninteracting case 

Let us now briefly address the noninteracting case, 
Ec = 0, where the Cooper pair sector decouples from 
the fermionic sector spanned by the |nin 2 n/) states. The 
Ec = 0 case with normal metallic leads instead of our 
dots was studied in Ref. [53] , while the dot case was also 
considered in Ref. |5^. An alternative way to reach the 
noninteracting limit is to let Ej —>■ oo, where all |iVc) 
states are strongly mixed and charge quantization effects 
are washed out [55]. In both limits, Ec —b 0 or Ej —)■ oo, 
one arrives at the same effective Hamiltonian H, where 
the parameters Ec, rig and Ej do not appear anymore. 

In the even-parity sector, iJ has a 4 x 4 matrix repre¬ 
sentation, which can easily be diagonalized to yield the 
ground state |4'). For the symmetric case with e = 0, see 
Eq. (31), and omitting the redundant Aj, indices, we find 


IvF) = aa{ (ey + ^e} + 8A2) [|000) + |110)] (47) 

- V8A[|101)-|011)]}, 

which is of the form in Eq . (|12| . We can therefore infer 
the concurrence from Eq. (]13[), 


turnover points between the CAR- and TP-dominated 


C(Aoo) = 


1 


\/r+^ 


Aoo = Vs A/e/. (48) 
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For e/ <C A, using Eq. (|^, we thus obtain the conven¬ 
tional exponential decay of the concurrence with increas¬ 
ing distance between the MBSs, C ^ ef /X ^ exp(—L/^). 
For 6/ = 0, the CAR contribution is precisely cancelled 
by the TP contribution corresponding to the second pair 
of states in Eq. (471, see also Eq. (131, and the concur- 


10 


rence vanishes identically. 

On the other hand, for e/ 3> A, we find (7 = 1 due 
to CAR processes, since entanglement is now dominated 
by the first pair of states in Eq. (47). The noninteracting 


system thus allows for CAR-mediated ideal entanglement 
as well, but (7 = 1 is reached only for rather small L 
where e/ » A, see Eq. For the realization of long- 
range entanglement, it is therefore essential to work with 
a floating island where Ec ^ 0. 

We conclude that the noninteracting case can also be 
understood quantitatively in terms of a competition be¬ 
tween TP and CAR processes. For even parity, we have 
seen that CAR processes dominate when a sizeable MBS 
hybridization e/ is present, which effectively corresponds 
to a rather short distance L between the MBSs. We note 
that for odd parity, the role of TP and CAR processes is 
effectively interchanged, again yielding (7 = 1 for e/ » A. 
Since we consider single-level dots and not metallic leads, 
otherwise competing processes such as normal reflection 
or local Andreev reflection El ED] are absent. As a con¬ 
sequence, the entanglement found in Eq. (48) is generally 
more robust. 


C. MBS hybridization effects 

For small L, the MBS hybridization e/ may exceed the 
charging energy Eq significantly. We then consider the 
parameter regime (23) but with arbitrary ratio efjEc, 


i.e., both e/ and Ec may be large compared to all other 
energy scales in the problem. 


Gate parameter rig = 1 


Let us start with the case rig = 1, assuming a symmet¬ 


ric setup, see Eq. (31). As shown before, for e/ ^ Ec, 
we then have (7 = 1 due to the TP mechanism. Now for 
finite (but large) £/ < Ec, the energy gap between the 
ground and excited state sectors in Eqs. (25) and (26) is 
given by AE = Ec — £/, see Fig. [^a). However, as long 
as the condition in Eq. (23) continues to hold with the 
replacement Ec —>■ AE, the situation remains conceptu¬ 
ally as before, and for ry = 0 (i.e., rig = 1) and Ec AE 
in Eqs. (27)-(30), we again find (7 = 1. 

On the other hand, for e/ Ec, our Ec = 0 discus¬ 
sion suggests that the CAR mechanism generates full en¬ 
tanglement of the dots with (7 = 1 as well. We therefore 
now consider a finite (large) charging energy with Ec < 
e/, where the previous ground and excited state sectors 
are effectively interchanged, Ho -o- Hi in Eqs. ( [2^ and 
(26). This situation is illustrated in Fig. [^b), where we 
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Figure 6. Color-scale plot of the concurrence for rig = 1 in 
the ef-Ec plane, taking Ej = A and a symmetric system with 
e = 0, cf. Eq. ([3l|). 


can read off the energy gap separating the two sectors, 
AE' = e/ — Ec- With the replacement Ec —>■ AE', 
we then assume that the condition in Eq. (23) holds 
again. In particular, the Josephson coupling Ej now di¬ 
rectly connects the respective states in the low-energy 
sector such that the CAR processes needed for generat¬ 
ing entanglement, cf. Sec. 1IIB[ are already captured by 
second-order perturbation theory, see Eq. Consid¬ 

ering the case ei = —62 and Ai = A 2 = A, and diagonal¬ 
izing the remaining 4x4 Hamiltonian in the basis (26), 
we obtain the ground state 

l^-) =aa{(f;j/2)[|ooo,o) + |iio,i)] (49) 

+ (r + Vr 2 + (Aj/2)2) [|110,0)+ 1000,1)]}, 


where the rate 


F = AVAE' = AV(e/-Ac) 


(50) 


connects jllO, 0) and jOOO, 1) through the virtual occupa¬ 
tion of excited states. Although the concurrence formula 
(13) is not applicable for |4t) in Eq. (49), (7 can be com¬ 


puted analytically, with the result 


C{Xf) = 


1 


1 + Xj 


Xf = 2T/Ej. 


(51) 


In accordance with Sec. [HTBl the concurrence vanishes 
for X/ 1, in particular for Ej 0. 

These analytical results are nicely recovered from a nu¬ 
merically exact calculation of the concurrence for rig = 1 . 
We show a color-scale plot of (7 as a function of e/ and 
Ec in Fig. confirming that lines of constant concur¬ 
rence are present when AE = —AE' = Ec — e/ is held 
constant. We took ei ^2 = 0 in Fig. 1^ such that both TP 
and CAR processes can generate st^l 
, In the crossover regime e/ 


le entanglement, see 
Ec, we find a (7 = 0 


Sec. HI 

line right at e/ = Ec, see Fig. This limit is correctly 
captured from the CAR point of view (even though per¬ 
turbation theory breaks down), since then the rate F in 
Eq. (50) becomes large and hence (7 —>■ 0 in Eq. (51). 
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Figure 7. Same as Fig. [^but for Ug = 0. 


Approaching this limit from the TP-dominated side is 
more subtle. The energy diagrams in Fig. |^a,b) reveal 
that there are six low-energy states at play for « Eq. 
These states are directly coupled via Ai _2 and Ej, which 
can be viewed as a mixing of TP and CAR processes. In 
that case, neither TP nor CAR are effective in generat¬ 
ing entanglement anymore since the previously excited 
states may now be occupied as well. 


1110, 0) are lowest in energy. We therefore obtain (7=1 
for almost the entire ef-Ec plane. The color-scale plot 
of the numerically exact solution in Fig. confirms this 
result, with the region AEq < A representing the only 
exception. In the latter region, one approaches the per¬ 
fectly disentangled limit e/ = Ec = 0, where CAR and 
TP processes interfere destructively. Finally, we remark 
that lines of constant concurrence are visible in Fig. 
when b" in Eq. (52) is held fixed. 


Concluding remarks 


Although the above results follow with the replace¬ 
ments Ec —>■ AE, AE', AEq in Eq. (23), we stress again 
that Ec itself has to remain large, Ec > |ei, 2 , Ai, 2 , ^^j|- 
Otherwise the description via charging energy parabolas 
as in Eig. is no longer sufficient and we expect devi¬ 
ations from the analytical results. Such deviations can 
be seen for Ec < A = Ej in Eigs. and where one 
effectively approaches the noninteracting limit Ec —>■ 0 . 


V. ENTANGLEMENT DYNAMICS 


Gate parameter Ug = Q 

Einally, we turn to the case rig = 0, where we consider 
optimal conditions for the CAR mechanism by choosing 
dot level energies with ei = —62 and symmetric tunnel 
couplings, Ai = A 2 = A. The corresponding energy dia¬ 
gram in Eig. [^c) shows two out of the six possible CAR 
processes, see also Sec. |IIIB| With the MBS hybridiza¬ 
tion e/, the hrst and the last state in Eq. (33) are still 
separated by the energy gap dEq from the ground state 
sector (32), while all other excited states in Eq. (33) have 
the gap^Eo = Ec + £/• Depending on the ratio eg jEc-, 
different sequences out of the six possible CAR processes 


then become important. Nonetheless, as long as Eq. (23) 
holds for Ec —>■ min^AEc, AEq), we may effectively ap¬ 


ply Eq. (35) again. We thereby find H as in Eq. (35) but 
with 


Up to this point, we have studied the stationary case, 
with time-independent parameters in iJ. In this section, 
we address the entanglement dynamics, (7(f), observed 
after a quench of a tunnel coupling where, say, A 2 is sud¬ 
denly switched on from zero to a hnite value at time 
f = 0. We assume that at times f > 0, we then have a 


symmetric system with e = 0 and A = 1 in Eq. (31), and 


we again consider the even-parity sector, V = 4-1. 

Studies of the entanglement dynamics can be very 
useful in revealing important timescales of the problem 
[m |49]. We will see below that similar to the static 
case in Sec.|IV[ perturbation theory is helpful in elucidat¬ 
ing the underlying physical mechanisms and the relevant 
timescales governing (7(f). In particular, our results for 
(7(f) can again be interpreted in terms of TP and CAR 
processes. 


h' = - 


1 


1 


Ec + e/ 2 ifc J 2{Ec + e/) ’ 


X^Ej 


(52) 


A. Finite charging energy: Ug = 1 


where we set rj = 0 (i.e., rig = 0). The first term 
comes from the two CAR processes with excursions to 
Q = ±.1 states, see Eig. [^c), while the second term de¬ 
scribes the remaining four CAR processes which also in¬ 
clude a virtual occupation oi Q = ±.2 states. Regard¬ 
less of which term dominates, throughout the regime 
AEq ^ max{X, Ej), we observe that CAR processes 
dominate the concurrence both for Ec 3> e/ and for 
e/ ^ Ec- This fact can be rationalized by noting that 
the inversion of ground and excited state sectors (dis¬ 
cussed for rig = 1 above) does not take place anymore, 
since for arbitrary ratio ej/Ec, the states | 000 , 0 ) and 


Let us begin with the case of finite Ec and arbitrary 
e/, keeping AE = Ec — e/ > 0, see Sec. IV C Eor sim¬ 
plicity, we put Ej = 0 and rig = 1, where we expect 
that entanglement is mediated by TP processes. We note 
in passing that in the presence of weak parity-breaking 
processes, e.g., due to quasiparticle poisoning, the even 
total parity {V = 4-1) sector is the true ground state 
at rig = 1. Suddenly switching on the tunnel coupling 
A 2 = 0 —)■ A at time f = 0 results in a time-dependence of 
the concurrence, where our numerical results for (7(f) are 
shown in Fig. [8 Here we start out from the initial state 
|(I'(f < 0 )) = 101 , 0 ), where the decoupled dot (QD 2 
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Figure 8. Time-dependence of the concurrence, C{t), ob¬ 
tained after a quench of the tunnel coupling, A 2 = 0 —>■ 1 at 
t = 0, for three different values of A_E = Ec — £/■ We con¬ 
sider the case of finite Ec with Ug = 1, taking Ej = ei ,2 = 0 
and Ai = 1. Exact results (black solid curves) follow from 
Eqs. (531 and (54l. The envelope functions (red dashed 
curves) in Eq. (551 accurately capture the slow part of the 
dynamics. 


in Fig. Q is empty before the quench. The post-quench 
state > 0)) is of the form 


|T'(t)) =ci(t)|110,0)+C2(t)|101,0) (53) 

+ C3(t)|011,0)+C4(t)|000,l), 


where the time-dependent coefhcients Cj(t) are subject 
to normalization, Cj{t) = 1, and obey the initial con¬ 


dition Cj (0) = Sj^ 2 - The solution is given by 


Ci(t) = -C4(t) = 


—iX 

V2n 


^—iEctl2 


sin(ili), 


(54) 


^-ieft/2 ^-iEct/2 


iAE 


sin(nt) + ncos(r2t) 


where = y/{AE/2Y -(- 2A^. The concurrence C{t) 
then follows by inserting the coefhcients in Eq. (54) into 


Eq. (13). Since for Ej = 0, total particle number is con¬ 


served, only the given four states in Eq. (53) with same 


total particle number as in the initial state are accessible. 
For that reason, Eqs. (1^ and (54) represent the exact 


solution for arbitrary AE > 0 as long as Ej = 0. Indeed, 
the Cit) results in Fig. are identical to those obtained 
from numerically solving the full Schrddinger equation. 

To gain physical intuition, it is instructive to compare 
these exact (but not very illuminating) results for C(t) to 
the corresponding predictions obtained from the reduced 
2x2 Hamiltonian H in Eq. (28), with the replacement 
Ec —>■ AE as in Sec. IV C| Using the parameters a and b 
in H, only a single oscillation frequency, ujp = 2\/a?- E 6^, 
appears within that approach, where for our symmetric 
system, we hnd ojp = 2X^/AE. The low-energy Hamil¬ 
tonian iJ captures only the slow part of the dynamics, 
i.e., the envelope function C'e(t). Moreover, its pertur¬ 
bative derivation requires the condition AE ^ A. On 
the other hand, for AE < A, we can directly separate 
the slow and fast variations in the full analytical expres¬ 
sion for C{t) described above, which gives Ce{t) also for 
AE < A. Combining both cases, we obtain a sinusoidal 
envelope function. 


C'e(t) = |sin(u;et)|, (55) 


with the frequency 


_ j AE/2, AE/X < 

“ I ^{AE/2Y + 2A2 - AF;/2, AE/X> 

(56) 

Notice that for AE/X < a/8/3, the slow part of the 
dynamics becomes universal (independent of the tunnel 
coupling A). For AE/X « y^8/3, C{t) is governed by 
two almost resonant frequencies and the envelope func¬ 
tion in Eq. (55) becomes less meaningful. Figure il¬ 
lustrates that away from this (quite narrow) transition 
regime, Eqs. (55) and (56) accurately describe the slow 
envelopes of all C{t) curves. For AE ^ A, the frequency 


We in Eq. (56) coincides with the perturbative scale Wp 


extracted from E[. In this limit, Eq. (55) becomes exact, 
C{t) —>• C'e(t), and the fast oscillations in C{t) disappear 
altogether, cf. Fig. By comparing the characteristic 
frequencies appearing in the exact solution for C{t) to 
perturbative estimates valid for AE 3> A, we thus have 
complemented our picture of the system in Fig. 
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B. Finite charging energy: rig = 0 


Such an analysis of the entanglement dynamics be¬ 
comes more difficult for other parameters within the 
large-ifc regime, We here briefly consider the case rig = 0 
with Ej > 0, where CAR processes can establish entan¬ 
glement. The Josephson coupling now makes the inclu¬ 
sion of higher-energy states with different total particle 
numbers necessary, and exact expressions for the enve¬ 
lope curves as in Eqs. (551 and (561 are more challenging 


to obtain. Nonetheless, analytical estimates are still pos¬ 
sible in the perturbative regime, where we obtain Ce(t) as 
in Eq. (551 but with Wg « Wp = 216" |, using 6" in Eq. (52). 


The CAR-induced entanglement dynamics at rig = 0 
is thus characterized by a sinusoidal envelope function 
again but has the rather long timescale ~ Eq/{X^Ej). 


C. Noninteracting case 


Another case of interest that allows for an exact solu¬ 
tion is given by Ec = 0. Since the Cooper pair sector 
decouples, see Sec. |IVB] we then have a set of four states 
describing the full dynamics again. The concurrence dy¬ 
namics after a tunnel coupling quench was recently ad¬ 
dressed for this limit of our setup in Ref. [25] , where oscil¬ 
lations in C{t) were reported as well. With e/ > 0, taking 
the pre-quench ground state |4>(t < 0)) = |il0), and ex¬ 
panding l'I'(t)) in terms of the four basis states available 
for even parity and Eq = 0, the post-quench [t > 0) state 
is as in Eq. (53) but with the initial condition Cj (0) = djp. 
By solving the time-dependent Schrodinger equation as 
detailed in Ref. [25|, C{t) follows again from Eq. (131. 

In fact, even though CAR processes are now responsi¬ 
ble for generating entanglement, finding C{t) for Ec = 0 
is fully equivalent to the TP-mediated case in Sec. |V A| 
after replacing AA = Ec — e/ —> e/. We thus obtain 
the same concurrence dynamics as in Eqs. (55) and (56) 
after this replacement (SD), including the results shown 
in Eig. |8| where now e/ « A separates the perturbative 


eTej 


regime (e/ A) from the regime where calculation of 
C{t) needs to account for all four states (cf < A), see 
Eqs. (|47|) and (|48|). 


VI. CONCLUSIONS 

In this paper, we have provided a detailed study of 
entanglement in a floating topological superconducting 
island hosting a pair of MBSs, each of which is tunnel- 
coupled to a single-level dot. The concurrence then pro¬ 
vides a convenient measure to quantify entanglement of 
the two dots. We have shown that this setup offers a ro¬ 
bust route towards the implementation of a “Majorana 
entanglement bridge” where the concurrence does not ex¬ 
hibit any decay with increasing separation of the MBSs. 
The only restriction for the separation of the quantum 
dots then is given by requiring a finite charging energy 


on the island, yielding much longer ranges of entangle¬ 
ment than in usual notions of exponential decay. This 
unique behavior is due to the intrinsic nonlocality of the 
fermion state representing the pair of MBSs. The under¬ 
lying mechanisms for entanglement generation have been 
identified in an intuitive yet quantitative manner. In par¬ 
ticular, we have shown that entanglement is created by 
the interplay of teleportation and crossed Andreev reflec¬ 
tion processes. With this understanding, we were then 
able to discuss the more complicated entanglement dy¬ 
namics after the quench of a tunnel coupling in terms of 
simple envelope functions. 

We have studied entanglement for the Majorana device 
in Eig. I^by employing the concurrence as entanglement 
measure. Although the concurrence has a rather formal 
definition, the predicted entanglement features should be 
observable in mesoscopic transport experiments, where 
the two dots are weakly coupled to additional leads as 
described, e.g., in Ref. m- Since current-current correla¬ 
tion functions can detect the violation of Bell inequalities 
[52] . signatures of entanglement are expected to appear 
in shot noise measurements, see also Refs. [53H55] for re¬ 
lated proposals. In principle also suitably designed con¬ 
ductance measurements could probe entanglement |56] . 

Recent experimental work indicates that 

devices corresponding to the model studied here are in 
close reach. For instance, charging energies of order 
Ec « 1 meV and TS gaps A « 0.2 meV have been 
reported in Ref. [15] , where quantum dots can be formed 
spontaneously near the TS wire ends when contacting 
them by Au leads m- In order to reach the regime 
A > Ec assumed in our model, the charging energy could 
be lowered by simply increasing the wire length. On the 
other hand, it may be possible to achieve significant en¬ 
tanglement between the dots even when A < Ec, where 
fermionic TS quasiparticles should be taken into account. 
When such states are localized in the bulk of the TS such 
that they do not have significant overlap with the Majo¬ 
rana bound states, they have no effect on entanglement. 
When low-lying quasiparticle states are localized near 
the TS nanowire ends, the physics discussed here may 
change, but entanglement may nonetheless be possible 
according to our preliminary analysis, see also Ref. [29j . 
Using NbTiN or InAs/Al [35] to proximity-induce su¬ 
perconductivity in a device similar to the one in Fig. 
extremely small poisoning rates have been reported as 
well. In any case, a detailed study of quasiparticle effects 
on entanglement is beyond the scope of the present work. 

To conclude, we are confident that future theoretical 
and experimental studies of similar setups utilizing such 
a Majorana entanglement bridge will benefit from the 
detailed physical understanding supplied here. 
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